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ABSTRACT 


A Carousel Remotely Piloted Vehicle (RPV) Recovery 
System consists of a vertical mast which supports a horizon- 
Du member that is free to rotate. The approaching RPV is 
Caught by vertical cables suspended from the horizontal 
member and the kinetic energy of the RPV is dissipated 
meroueni the motion of the recovery system. This thesis 
presents а simolified dynamic analysis to describe the 
motion of en REV after impact witn a recovery system simu- 
Having the carousel system. The equations of motion for 
the RPV were obtained from Lagrange's equations and modeled 
using the Continuous System Modeling Program on the IBM 
360-57 computer at the W.R. Church Computer Center, NPS, 
fhe results showed that damping was required on motion of 
Pome спе horizontal member опа the suspended cables in order 
GO prevent possible damage to the 150 pound EPV at an as- 
БИЛЕ 50 knot approach speed. А sensitivity analysis was 
performed ру varying the system desizgm parameters and the 


Mesulos are presented in the body of the thesis. 
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па ЕРА РОВІ 


Extensive efforts on Remotely Piloted Vehicles (RPV) 
have been underway for many years, covering a wide range of 
Soolicatlons and missions, Ref. 1. It is only recently, 
however, that both the need and the technological capability 
for mini-RPVs have coalesced to the point that the develop- 
ment of these systems has become imperative. A mini-RPV 
is generally defined as an RPV that veigns less than 200 
lbs mission loaded. 

Àn essential part of an RPV system is the recovery 
system. The Army has recently experienced several mini- 
RPV losses ave eres аатестту ог indirectly to The recovery 
method. The methods used to recover mini-RPV's vary from 
catching the RPV with an arresting cable similar in prin- 
шоље to that used on an aircrait carrier, to guiding tne 
плита ПРУ сілтестіу into a vertical net which is designea to 
flex and absorb the kinetic energy of the RPV. A completely 
setisfactory recovery system has not yet been established, 
hence studies into alternate recovery systems are of con- 
tinuing interest to the Атту. 

О Шр се Os wune analyses described hercin was to 
model the dynamics of one alternate recovery system candı- 
date, the "carousel" method. The carousel method was 
Mao сесси ас a candidate for є recovery procedure 


in a concert study performed by Teledyne Ryan for the Army, 


Ti 








ШІ 2 pub the initial concept formulation did not include 
significant dynamic considerations. Therefore, the analyses 
Eouducbed іп this thesis, which are primarily oriented to- 
wards system dynamic response behavior, should serve in the 
ог providing additional information ior concept feasi- 
bility considerations. 

The carousel recovery system, shown in Figure 1, was 
concevtualized in this analysis as being composed of a 
vertical mast that suvported a horizontal member which was 
allowed to rotate about the vertical mast. Hanging from 
the horizontal member are cables that will be caught by 
the wing tips of the approaching RPV. The RPV's kinetic 
energy will be dissipated by the recovery system through 
ШЕЛІ ола means that will be evaluated in this analysis. 

The REV must be brought to rest without allowing 
excessive ВРУ swing since excessive swinging motion could 
result in the RPV hitting the recovery system structure. 
This consideration was of paramount concern when evaluating 
mme merits oi recovery system parameter variations. 

Por the purpose of this analysis tne mini-RPV to be 
recovered was assumed to be the AQUILA, which is a 0.5. 
Army mini-RPV developed by Lockheed Missiles and Space 
пазл, Rel. L, end 1s sketched in Figure 2. It has the 
molto mias nominal characteristics: wins span, 12 ft; 
weight 150 10; approach speed 50 knots; maximum loading, 


ес о edel, 32 lateral. 


152 
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А. ASSUMPTIONS 

simplifying assumptions, as defined below, were made 
when establishing the equations or motion using the pro- 
cedures suggested by Lagrange's equations, Ref. 3, in order 
ОШ оеп the basic principles of the carousel recovery 
system dynamics. Inclusion of system details such as re- 
covery boom elasticity and RPV mass moments of inertia were 
considered as items which might tend to mask the basic 
question of feasibility at the present stage of analysis. 
Under these assumptions, the results supported the viability 
of the carousel recovery method and a basis was established 
upon which further analyses may be performed when the sim- 
peer vins assumptions are removed. 

The recovery system was initially assumed to be con- 
servative, that is, friction losses were ignored and as 
soon as the RPV contacted the recovery system the engine 
was secured. This assumption was subsequently changed by 
the addition of damping to the recovery system configuration. 
The horizontal member was assumed to be rigid. A careful 
choice of material and structural characteristics for the 


horizontal member can minimize its dynamics and allow it to 


be aporoximeted oy a rigid member. The cables hanging from 


і 


the horizontal member were replaced by massless thir 


covery rods. The RPV was treated as a point mess and was 


ty 
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аво DY one recovery rod. Treating the RPV as a point 
mass simplified the coordinate system used in the analysis 
by allowing the mass moment of inertia of the RPV about its 
pusncspal axes to be neglected. The RPV was also assumed to 
approach perpendicular to the plane formed by the horizon- 
tal member and the vertical recovery road that caught it. 
The approach velocity of the RPV was assumed to transfer 
undiminished to the recovery rod giving tne rod an initial 
ensular velocity only in the D direction. Reference should 
be made to Figure 3 where the orientetion of the generalized 
coordinates ©, D and и is portrayed. These coordinates 
are a form of Euler angles, indevendent of each other, that 
together completely specify the position of the RPV relative 
to the fixed reference frame x-y-z. First @ defines the 
rotation of the moving coordinate system x'-y'-z' about the 
fixed axis 2. It describes the rotation of the horizontal 
member in the x-y plane. Next (p defines the motion of the 
ШОШО ОП x=” within the moving coordinete system. This is 
portrayed in Figure 3 as the motion of the RPV projected on 
ШООК у'—2' plane. Finally, W depicts the remaining motion 
of the RPV after it rotates through Ө ала Ф. 

The system identified in Figure 3 is the reference or 
base configuration that was assumed for the anelysis based 
on judgement of the RPV's weight and size. In this refer- 


ence configuration the horizontal member, R,, was 28 ft 


I 
long and made of 6-inch aluminum I-beam having а mass 


moment of inertia of 1159 ft_-1b-sec“ АОС Рао ОУ 


16 








Figure 3 Simplified Recovery System 
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Une recovery rod R., was 6 ft long. The RPV of mass 4.66 


p 
slugs (150 lbs) was suspended at the lower end of Ry. With 
an assumed approach speed of 50 knots (84.4 ft/sec), the 

initial angular velocity, Фо of the massless rod Но was 


14 rad/sec according to the relation of equation (1): 
. Approach velocity (ft/sec) 
Do = Ro (E == — (1) 


A linear, torsionel soring, of spring constant K, lo- 
cated at "0", the pivot of the horizontal member, was 
initially used to resist the rotation of the horizontal 
member. without any damping in this reference configura- 
tion it was further assumed that the spring's stored energy 
could be eliminated from the system through a ratchet re- 
lease mechanism, after the horizontal member had swung 


through its initial maximum displacement. The ratchet 


mechanism was not modeled, as it was felt, and later con- 
tirmea oy the results, that the critical motion of the RPV 
would occur shortly after impact with the recovery system. 
Simulation of the ratchet mechanism would have no effect 


on these results. 


В. OVERVIEW OF ANALYSIS 

The general approach of the analysis was to model tne 
RPV motion after engagement with the recovery system using 
the generalized coordinate system Ө-ф- 1 indicated in 
Pig. 3. ii the assumptions orf а rigid horizontal member 


and a point mass RPV representation were removed, a minimum 
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of five more generalized coordinates would have been re- 
фе ес (је equetions ot movlon for the Three degree 01 
freedom system were obtained using Lagrange's equations 
tor a conservative system. From these equations the 
acceleration, velocity and position of the RPV as a function 
of time were obtained in component form along each coordi- 
nate direction through the use of the Continuous System 
Modeling Program (CSMP). The input parameters were then 
varied and their effect on resulting dynamic motion studied. 
Changes to the basic recovery system which included the 
addition oi damping terms were later entertained as a con- 


-Еошедсе от Initial results. 


С. DETAILED ANALYSIS | 
XL 5u-cranse'svguavLon 
The equations of motion were established using 
i шешее'5 equation, Ref. 3. Three equations oí motion 
were obtained Tor the system, since there were by assump- 
он only three generalized coordinates: the coordinates 
а O MEN en ee 
O, D and Й, аз defined in Fig. 3. Lagrange's equation 


mame Conservative form, without any generalized forces, is: 


af oat дт ð 


= E —— + =m M — Os gps == ЯМ 2 22) 
at e ce El у 4 ит nd 
902 = "xl Ca 
anere 
IE = г ali | 1a zt = ma == 
a generalized coordinate (qi 6.2, ORLE Y) 
Qi = time derivative of q, coordinate 


ЩО 





= 
|| 


H = Kinetic energy (ft-lb) 


=, 
|| 


V (az) = Potential energy (ft-lb) 


The kinetic energy was defined relative to the inertial 


reference frame by the following expression: 


T = BI 0 Enny (3) 
where 
io = Moment of Inertia of the horizontal member 
about pivot point "0" (ft-lb=sec”) 
il = Mass of RPV (slugs) 
Ө = Angular velocity in O coordinate direction 
(rad/sec) 
Увру =% Velocity expression for RPV derived in the 


певен о оно ван Sa) 


The potential energy included the gravitational term acting 
upon the RPV mass and the stored energy in the torsional 
constraint spring. Equation 4 expresses the potential 


energy relative to the inertial frame of reference as 


ош о о: 
V= МеВ, (1-cosoos P)+ix 6° (4) 
where 
Ro = length of recovery rod corresponding to RPV 
distance from the horizontal member following 
recovery, (ft) 
Е КОС а of spring at pivot of horizontal 


member (ft-1b) 
20 





Face its exactness is critical to the evaluation of 
Lagrange's equations of motion, the RPV velocity is treated 
separavely in the ensuing section using the generalized 


Booardinates. 


Pea Derivation of REV Velocity 


ООО ЖЖ ОЕ pl ss Р hs PSD 


То derive an expression for the RPV's velocity 
EBbPer impact with the recovery system, the position of the 
RPV within the moving x'-y'-z' coordinate system as & func- 
tion of the generalized coordinates 9, Wwas transformed 
to the fixed inertial reference system x-y-z. The time 
Кави ваме ош Loe RPV'S position, expressed as a function 
of the three generalized coordinates Q,@ and ү and ref- 
erenced to the fixed inertial system, represented the ab- 
-solute' velocity of the RPV. 

The RPV's position within the moving frame x'-y'-z' 


was expressed as follows: 


ХІ = -R, sinW (52) 
y! = Ro sinfcos Y (56) 
Ди 595 совфсоз у (5c) 


МЕ пас оте о to che теа inertial coordinate 


system may be stated as: 


x cosQ - sing -R sin Y 8. со50 
= + 
y sinO  cosQ R,sindcosy Boma 
- + (ба) 
= t me a a 
Z = Z = -R,cos@cosy (въ) 


Es 





In an alternate form, the three components expressing the 
absolute position of the RPV mass following recovery en- 


gagement may be expressed as: 


х = В созб -R,sinWcos® -R,sin® сов /віһб (Ta) 
y = R4cosQ -R,sinYsind +R sing cos Ycos 8 (7b) 
2 = -R,cos@cos ү (ЛО) 


The absolute velocity of the RPV мав obtained by taking tne 


first time derivative of equation 7, thereby yielding: 

x = [-R sing +R, (sin Wsin 0 -sinpcos Wcosf ) Ө 
-(R,cos® cos YsinO Io +R,(sin® sinWsin@g- 
cosWcos@ ) Y (8a) 

y =" | В соз6 -В (в1и YV'cosA +sind)cosVsinA) 

қ р Ө 0)6 
+ (В созфсоз /совб )ф -R, (cos ern RSL ene 
соз 6 (80) 


Ro (віпфсовзұ/)ф +(соѕф ѕіп ү) y] (8с) 


С. 


ое Ц In Lasranzets 
2 б 
eguation, Увру , was obtained by squaring the components 


Ss marasely then addins them; 1.¢., 


2 PAT 


ө > . 2 = ща S 
Mo expressions гот x, y and z are equations 10а, р 


End c. 
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|-R sing +R,(sinYsin Q-sinbeosWeosQ )] 2 2 

+82с05° феоз? зат" 0 ф“ 

«Ro (sin sinV'sing -cos V/cos Q ) V? 

-2R,cospcosY/sing [-2, sin g+R,(sinysing -sing 

cosyrcos8 ] 6 Ф 

+282, (singpsinysinG -cos/cosg ) [-2, sing +R, (sin 

sing -sindcosyreosQ ) 0 y 

-2R, соефсозу віп 8 (зіпф віп У віл 0 -созу сов Joj 
(102) 

R cos -R,(sinYcosg +sindeosWsing y] OHR,“ 

cos фсов г сос“ o" 

27" [eosyrsingrsinfsinycosg] * y^ 


+2R,cos Фсозу/сов 9 (2,2950 -R, (sin Wcos Өззіп 
cosWsin® Je D 


-2R, [cosYrsing +sin(Dsin cos 9 [Ry cos -R, (sin W 
cosQ + 511 фсов Ива 9) Ө 

-28> созфосоз Y cos - [cos Wsing +sin@siny 
cos 8| у D 


(10b) 


ГӘ 
{АЈ 





22 = ES bin p cos y p +сов ф sin” yr + 2в1афсовф 
sin сос Фф (196) 


aD © 
Шеп х,у”, 22 were combined, several cancellations 


| 2 < T 
occurred and the final expression for V ео: 


REN 
2 Ze EE 2 В а 
У Ру = [Ra +R, Sin V+R, sin (D cos Y -2R,Rosiny | 9 
20 Ne г „Ж 
+8, со уф +В. Y +R, cosp cos Y(R -R siny) Oo 
а < Ви T 
+(-2R,R sing siny+2r, “sin )0 y 
(11) 
кри Бег ла об о: Equations of Motion 
а. б Equation of Motion 
Now that an expression тот Em has been 
obtained, Lagrange's eauation can be evaluated according to 
section 11.0.1. The equation was evaluated for each gener- 
alized coordinate Ө, HW and. The O eauation of motion 
was obtained оу evaluating the Lagrange's equation using 


the @ generalized coordinate. 


4% |96 9 20 = 
Asain: 
aa я . T 2 йрт 2 
Kinetic energy, T=5I, 0“+3lVppy Gis) 
ee 2 
Potential energy, V—lisR 1-созф соз +-К © 
1 | (14) 
97 
22-0 (е) 
20 
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г - 2 АА DD 2 А 
42%) -Т 0 +М KR, +R, sin yrtR, sin (cos y -2R,R,siny)0 
+Rocos@cos V^ (R4 -85 sin Wd +(-RyRosinQ siny+ 


Ro sinp) Y 

+(2R, sinyreos Y y +2R, sing cosp cos y d 
-2R,?sin?p eos sin умов, В,созуг ур б 
-R,sinDcosY(R, -R,siny) H -R,sinyrcos PIR, -R, 
зе Фу 

Ro cos усозфоозу Vr -RitecosQ sind y-R122 
sinfcosY y. 

+В созф oy (17) 


Неа та шале ecuations 15, 16 and 17 into equation 12 and 
(КЕЛП тп through by (ив г) yields the following Geguation 


SIN motion: 





I R » 
2e (gh-sin V?) сов віп'Ф )6 
TR 2 
2 





R n" 
квіпф(1-алвіпд) уг 
M 
в. : 
+cosp (s=cos//-žsin2 Y) D 
2 
2 На ° А 
+sin2W cos D -2cos) Ye 
ме 
+5112ф cos Yr ф 6 
ty 
+cos@ (- -aysin ny -cos2Yy +1) p Фу 


В 
-віпф (qecosi—ssin2y) ф і 


Ry 
= — 51 n@cos Vr V 


a 


Equation of Motion 
Similarly the (D equation of motion was ob- 


tained trom the following Lagrange's equation: 


1 д 
3 (8)- ді + 9V — 0 (19) 


where 
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-2R cosy siny Wo-R,sindeos YiR--R,sin 1) Фе 
-R,sin /созф (®{-Е„зїпүг)Өүг -В„°созф сове ye) 


(20) 


E = М [a зіпфсозф cosyg"-2,singp cos (В -В siny) O 

+(-R1Rocos&® siny+Ro“cosQ) )6 y] (21) 
IV - MR singcos (22) 
ID gu 


Substituting equations 20, 21 and 22 into equation 19 and 
EN... — 2? 2 
dividing through by UR, ) yields the following (D equation 


of motion: 
На 3; 
g-cosqQ сову -cosD sin2Y| 6 
2 
2 .. 
£cos^ У Ф 
-sin2Y YG 
-созф (соз2\/+1) Оли 
-isin20 cosy д? 


т sin coswe= 0 
Eo Ф E (23) 
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S yr Equation of Motion 
ade W equation of motion was Similarly obtained 


using the following Lagrange's equation: 


21 t QC 
ea) E c 


where 


а [Эт DEM іх о, | ске! ee 
== E) = R, M (= псозф siny+cos@) D6 - pon" 
R 


СЕ ni siny+sing)0| (25) 


A R, “и (siny cosy -sinp cosysiny - тоу) 8^ 
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Substituting equations 25, 26 and 27 into equation 24 and 


Anding through by the (HR,“) yields the following W 


equation of motion: 
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4. Solving tne iquations of Motion 
The equations of motion were solved for time his- 
tory information oi the RPV's position, velocity and accel- 
eration using the ТВЫ Continuous System Modeling Program 
(CSIP) which was available as resident software on the IBM 


360-67 computer located in the W.R. Church Computer Center, 


Naval Postgraduate School. The use of CSIP as a high-level 
programming tool is quite similar to the familiar engineer- 
ing language of FORTRAN, ano programming details may be found 
Beine test by speckhart and Green, Ref. 5. Details of tae 
Men program usea are presented in the Appendix. To use 
this software package, the three equations of motion were 
first solved for the highest order derivative Ө, Ф and Y 
її terms oí the three generalized coordinates and their first 
time derivatives. Then, in a manner similar to an analog 


computer solution of differential equations, tne COMP program 


proceeded to integrate these equations numerically for values 
for 0, D, W, 0,  snà Y, which were fed back into the 


equations for 0, H ana W. This iterative method ої 
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solution was repeated at each time increment. The resulting 
values for 0, Ф апа У were listed as a function of time 


and from these the motion the RPV was studied to assess the 
feasibility of the recovery system when a particular set 
of structural parameters and initial RPV conditions wes 
Specified. 

Since the governing equations of motion (egns. 18, 
23 and 28) were inertially cross coupled, it was necessary 
to rearrange the equations into an inertially uncoupled form 
that would be accepteble for solution by the CSMP. In matrix 


won, equations 16, 23 and 28 may be expressed as; 
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substituting equations 30a through 30h in matrix expression 


29, 9, Ф and W were solved using Cramers' rule, Ref. 4. 


mE | БОД, 
р IA (31a) ЈЕ (310) 5-13 (310) 
where 
ШЕ 12) 758 
Я = |221 9 893) = 812 423 831-311 823%821 213 
| 231 аз2 9 (32а) 
p TERRE MESES 
Al = |» О 525702 21: 
Е O T 





Е 
[42| = | 221 Po 223 |= 231(91273-b281 3)-b3(8,,293-2), 31 3) 
par * PU (32е) 
aL ts E 
EE ау 0 92 |- аҙ1(92212)-(92311-04221)%9 
| a3] 1 m | (-821242) TT 


Squations 30a through 32d were nrogrammed on the IB 360-67 
Euer usina СоОМР, «пе details of which are described in 
tne Appendix. Tne program evaluated 6 ; d У у , 6 ; ф : Y ; 
END. and yo as a function of time from the initial con- 
dition of RPV impact. The evaluation of the computer 
generated time history results is described in the next 


ото. 
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ТТТ. SYSTEM DESIGN ANALYSIS AND RESULTS 


ІІІ UNDAMIPED SYSTRH 

The reference or base configuration of the recovery 
eren described in section II.A. and Fig. 3 was the first 
configuration considered in the analysis. From this con- 
mesuravlon each parameter was varied individually for various 
RPV apvroach speeds to assess the system's sensitivity to 
their changes. 

From the outset of the investigation it was felt, as 
ЕШ Шектес by carly results, that the motion ої tne RPV 
along the (D coordinate would be a primary system constraint. 
imourdeline value for D maximum (Dmax) was chosen to be 90 
degrees, which corresponds to the recovery rod (Ro) reaching 
the horizontal. For the recovery system to pass its first 
eoa renent, an REV approaching at 50 knots would have to 
be recovered with the maximum amplitude of D not exceeding 
90 degrees. 

The results of modeling the reference configuration 
suggested that energy disipation would have to be included 
ШОКОТ лт Qrar dorn to 90 degrees or less. Without any 
damping the system exhibited unacceptably hign magnitudes 
TOT Ф that corresponded to the RPV completely rotating 
about the horizontal member wnen approach speeds exceeded 
Б 055. THis can be seen by referring to Table I where 


ЫС sted tor whe Dase configuration as a function of 





TABLE I 
VARIATION FOR UNDAMPED SYSTEM 
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Mime alter impact for avproaen syeeds rengins from Y knots 

( Ф,-2 паси сес ner, Equation 1) to 50 knots ( D,=14 rad/sec). 
more that Tor db, =1 ео. or approximately a 25-knot 
approach sveed, tne magnitude of @ at 0.7 seconds after 
impact was greater than 180 degrees. The REV has started 


ootate over the horizontal member. Variation of the 


System parameters did not correct tnis deficiency. 


Б СА ТРЕНС ADDED ТО THE HORIZONTAL MEMBER 
Поа ел ок го correct tae excessive RPV swing in the @ 


coordinate direction, a Rayleish dissipative energy term, 


D, Was introduced into Lagrange's equation: 


ен с - о 05, 
994481 Ci 90i 944 (33) 


2 
where D=3C,0 


This additional term precipitated through the equation 
ӨШ motion derivation of section II.C. with the addition of 
emily a 0,0 term to the 6 equation of motion, equation 18. 
(hts enanze vas easily incorporated into tne computer pro- 
Eu ues indicated in the Appendix. 

An initial estimate of the magnituce of the damping 


constant C4 Во и де оке општој тес Linearized 


penculum eauetion pelov. The damping value wes varied to 


— 


cetect the system's sensitivity to "E 
6+2 2w.6+@, 0= 0 
о еы (34) 
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A (беу = а 2 fa 1 p 
Por values of 1,=1159 1b-ft-sec , K=50 £i-lo and б- 0.15, 


Мол пас selected to introduce just a small amount of 
damping, the initial estimate oi C4 was о 


The results of adding damping C, to the pivot of the 


ab 
horr2onuval member in the base configuration indicated that 
it Was only partially successful in reducing Qmex. This 
can be seen by referring to Table II where @ is listed for 
MECON guration as e function of time efter impact for 
“Шр бос speeds ranging from Y knots, where ф,-2 rad Sec. 
to 50 mc v aere &,=14 red/sec. Note that the magnitude 
oi D was still avproaching 180 degrees for approach speed 
Pomel conditions in excess of 35 knots ( &,=20 rad/sec). 
Variation of system parameters including C4 Uoc De 
то turtner reduce Фпах, as shown in Figures 4 through 9, 
упеге ihe base configuration of 24220 ft, H5-26 Ft, 11=4.66 
slugs, 1,=1159 ft-lbs-=sec", K=50 ft-lb and 0,=75 ft-Lb-sec 
Мт осшес win only the indicated parameter varied, ceteris 


paribus. 


С, DAUPING ADDED TO RECOVERY ROD 
In a further effort to reduce ine maximum swing, Фтах, 


of the recovery system, damving was also adcec to the 





TABLE II 
VARIATION WITH Cy DAMPING 
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recovery rod. The added camping was included in the system's 
mocel, without cross coupling, oy modifying the Rayleigh 
dissipative energy function to the form of: 


ие 
D=3C,0 “+350, D (35) 


The additional dissipation term precipitated through the 
equation of motion derivation described in section II.C. as 
the addition of a 0,0 term to the @ equation of motion, 
Pellauilon 23. This change was also easily incorporated into 
the computer program as indicated in the Appendix. 

An estimate of Со was obtained in a manner similar to 
that described in the previous section using the simple 


pendulum equation: 
.. ч 2 с | 
Фф 2СбоуФ% шаФ- 6 (36) 


where 


Only a slight amount of damping, Co=3.24 ft-lb-sec, 
added to the recovery rod of the base configuration, in 
addition to tne horizontal member damping, was sufficient 
to reduce (max to approximately 90 degrees at RPV approach 
speeas of 50 knots ( Ф,-14 rad/sec). This can be seen by 


referring to Table III where @ is listed as a function of 
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time after imvact for approach Speeds ranging from 7 knots 
( &,=2 Rad sec MS ( ф,-14 rad/sec). The maximum 
magnitude of б апа и for the modified configuration at the 
50-knot RPV approach speed was 50 degrees end 63 degrees, 
mes pectlively. 

Further study of the modified base configuration vas 
Mer rrormea to obtain the magnitude of the RPV's linear 
acceleration. The absolute value of the RFV's linear 


acceleration was obtained using the following equation: 
в РЕТ 
Acceleration =y X"+Y"*Z Cay) 


woere x, у апд 2 аге defined by equations 7a, b andc. A 
plot of equetion 37 as a function of time for tne modified 
Baer zeonfiguration is presented in Figure 9. It may be 
noted that the masnituae reduces to below one "g" within 
one second aiter the RPV contacts the recovery system. 

Tae sensitivity of (max to each of the system para- 
meters is indicated in Figures 10 through 16, where the 
Pease contiguration of R,=20 кре R,=6 Tt, 11=4.66 slugs, 


2 


METIDO Ti=lb=sec", K=50 ft-lb, C,=75 ft-lo-sec and C 


> 
3.24 fi-lb-sec is plotted with only the indicated parameter 
Varied, ceteris paribus. It can be noted from Figure 12 


that the system was sensitive to the recovery rod length, 


R which may pe interpreted as an implication tnat the 


2? 
ee e ement point for REV capture on the vertical cables 


fas considerable importance to system dynamic response. 
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The system was also sensitive to the value of SE but 
more sensitive to Со: This, however, should not present a 


mrobilem, as 54 and C, will be set by the design and will 


2 
not vary in operation. 


The system was relatively insensitive to variations of 


$ 


the point along the norizontel member, Rj, where the RPV 


was caught. The system was also relatively insensitive to 


varlations in the moment of inertia, То, өш the horizontal 
member and the mass, fl, or the RPV. The maximum displace- 
ment of Ф was not affecte at all by changes in the spring 


constant, K. 


D. IMPLEMENTING DAMPING ON A RECOVERY CABLE 

The damping desired at tne pivot of the horizontal 
member can be readily implemented TmsaeLlalzeonescruervıon, 
Ше ег, 1% was not readily anparent how the damping of 
mie vertical rod, which was a cable in the actual system, 
some бе implemented. 

Two Simple pendulum designs vere analysed and modeled 
using the same techniques described to study the recovery 
system. The intent was to see if pivotal damping could 
EE теч effectively to control the swing of a cable. 

The tirst design is shown in Figure 17 and simulates 
a cable design for the vertical massless recovery rod that 


was used in the recovery system analysis. With no pivotal 


damoing, an initial velocity was avolied to “2 to mak 


M approach a high angle without exceeding 180 degrees. 
1 4 


24 





Cable 


М2 


Figure 17 





Pendulum Design A 
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Then Cemping was epplied to the pivot of the cable support 
M, to judge its effect on reducing the magnitude of Ф,. 
The results are plotted in Figure 18. This approach of 
aoplying pivotal damping did not efrectively reduce oe 
pomcesired. 

Another design, shown in Figure 19, was also analysed 
morsee if it could, through pivotal damping, control the 
maximum magnitude of Po: Again the design simulated a 
cable instead of tne massless recovery rod used in the 


analysis. A spring was used at the pivot to give M a 


| 
near horizontal static deflection. With no damping applied, 


an initial velocity was found for My that caused Domax to 
approach 180 degrees. Pivotal damping was then applied 
to MH, to observe its effect on the magnitude ог Фо. The 
results shown in Figure 20 indicate that pivotal damping 
applied to this pendulum design was also ineffective in 
reducing Open. 

While neither of these designs demonstrated effective 
Sontrol of Do» tney by no means exhaust the possible 
metnods of implementing damping in a cable design, and 


further study is possible. 
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Figure 19 Pendulum Design B 
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IV. CONCLUSIONS AND RECOMIENDATIONS 


For the carousel recovery system to be feasible at 
approach speeds of 50 knots, a relatively small amount of 
damping was required at the pivot of the recovery rod in 
addition to damping at the pivot of the horizontal member. 
ihe damping at both vivots was reauired to prevent the 
maximum swing of the point mass RPV from exceeding approx- 
imately 90 degrees in the (PD direction, the coordinate that 
exnibited the largest change in magnitude. 

the recovery system was sensitive to the vertical 
height of the RPV at the time it was captured. The greater 
une distance, Eo) between the RPV and the horizontal member, 
the greater the swing or the RPV relative to the horizontal 
member. when selecting parameter values such as the damp- 
ing constants under the most severe conditions, tne RPV 
snould be assumed to have contacted the recovery system at 
the maximum length of the recovery rod or caoles. 

The recovery system was relatively insensitive to 
Variation in the distance the RPV was from the vertical 
mast wnen it was caught. The system was also relatively 
insensitive io variations in the mass of the RPV and the 
moment of inertia of the horizontal member. 

Variations in ihe spring constant for the torsional 
Soring located at the pivot of the horizontal member had 


ша и шелеттесъ om tne Swine ot the RPV. Since it was not 
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mayortanu to controlling the motion of the RPV, it could be 


removed from the design. 


The masnitude of linear acceleration exhibited by the 
RPV, with the recovery system design that incorporated both 
horizontal member and recovery rod damping, was within 
Aquila load limits. 

Based on this enalysis the carousel recovery system is 
still viable; however, it is apparent that tne system will 
be Sensitive to approach techniques used during the recovery. 

It is recommended thet this analysis be continued, in 
В ет to treat the RPV as a distributed mass with its om 
moment of inertia and coordinate system, and caught by two 
ШЕ Сту rods or cables, опе тог each wing tip, instead oí 
one. Expanding the scope of the recovery system model would 
Gao нае а more detailed analysis of the RPV's linear lateral, 
longitudinal and normal accelerations ana hence insure that 
Ley were individually witnin the RPV's design limits. 


ft is also recommended that the effects of an RPV 


y 


2200roaching at angles to tne recovery system be investi- 


С 


le 


Paved; this should result in increased motion in the Y 


C 


direction. 
It is further recommended that continued study snould 
БЕ спе ol a means oí avoolyine dempine to the swine of tne 
DT O D => сэ 


Метал сал recovery cables. 
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АРЕ ОТА 


MODELING THE EQUATIONS OF МОТТОМ USING CSHP 


~“ 


Ап application-oriented software package called the 
Continuous System Modeling Program (CSMP) was used to model 
and solve the nonlinear, cross-coupled eauations of motion. 
CSIP is a powerful software package that can handle this 
Be oT problem routinely. In the analysis Der Ommed sin tnis 
thesis two integrations using a fourth order variable inte- 
gration step Runge-Kutta technique on each of the three 
Баша ойе от motion, Her. 5. This sophisticated integration 
technique has the advantage of automatically adjusting the 
time increment of ЫЕ шл ко лее кые demands of the 
ШШепіс conditions oi the Simulation. The absolute value of 
the estimated integration error and the relative magnitude 
or the estimated error are compared with user specified 
error—bounds. The step size is then adjusted to meet the 
desired error criteria. In the analysis performed in this 
presis tie default value of .0001 for both the absolute and 
ete ive errors was used. 

The program structure ol COMP is composed of three 


segments; inivial, Dynamic, and Terminal. Generally, data 
statements will appear in the Initial ас эш сг 0° 206503; 
сешосштастотв Слаб are required to ре performed only one 


time during tne simulation can be conveniently placed in 





this segment. The Dynamic segment is usually composed of 
structural statements that describe the dynamics of the 
Een or explicitly deseriíibe a set of differential equa- 
tions. The Terminal segment is the last segment in the 
program and is usually made up of output control statements. 
The actual program used in this analysis for the un- 


damped system will now be described. List of symbols used 


ишле Computer Programs! 


RL - Length or horizontal member in feet 

R2 - Length of recovery rod in feet 

М - Mass of RPV in slugs 

MIO = Moment of inertia of horizontal member in ft-lb- 


sec“ about its pivot point 

K - Spring constant, in ft-lb/red,.of spring at pivot 
ОТ ПОТ опса шере 

ENT - initial angular velocity of the recovery rod in 
the @ direction in units of rad/sec also equal 
to approach velocity of RPV in ft/sec divided 
by Re 


G ссе еа отоп оиса зо ft/sec! 
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XADOT — 
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DOT 


у 
C 


X6DOT - 

CO - 64: damping constant in It-lb-sec for the damper 
located at the pivot of the horizontal member 

С02 - C,, damping constant in ft-lb-sec for the damper 
located where the recovery rod attaches to the 
horizontal member 

AD AN ever nan гот ее петле in Section ir.0.4 

ADI - Determinant of Al as defined in Section II.C.4 

AD2 = Determinant of A2 as defined in Section ТТ.С.4 

AD 3 сонне та о ено = се пе а ест от. Е СА 

BI - is b, defined in Section I1.C.4 

B2 - is by defined in Section I1.C.4 

33 - is by defined 39899060109 17.9.4 

Al - is 241 сего соо 1.0.4 

А2 - is a,, defined in Section II.C.4 

АЗ - is 24 3 del ino сео слог 11.6.4 

AL - is 


253 erica ec Lon I1.0,4 
A; B; 2 D, Ша n He ШЕ 5 T = were used DO Simpliry 
program notation and are 


defined ір the program 


Be soription or. ıne mrosran shorm in Figure 21, follows: 

Lines 1-3: Mac a ша ве стене, parameters 
ешо eelse mtemee here. 

Lines 4-34: Maze Но апе Dymamic segment. 


kines 2-15: иис Oba 1011. 
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Figure 21 Computer Program 





Lines 14-28: The algebra used to solve for 6, Ф and Y 
Prior 50 mam one as mCeSseri pe: in 
DOCU lO a 

Lines 29-31: Perform the first integration from angular 
accelerations to angular velocities. 

Lines 32-34: Perform the second integration from angular 
velocities to angular position. 


Lines 35-41: Make up the Terminal segment. 


Mime 36: Limits the time period that will be modeled. 


Lines 30-39: Prints output as a function of time. 


When damping was added to the horizontal member the 


only changes that occurred to the program were the following: 
Line 2: за ООК 
Line 23: add  -CO * XA. 
When damping was also added to the recovery rod the 
additional changes that occurred to the program were as 
MOS отс: 
Mne 2: зоо CNN 


ШІЛе6е 24: add  —-CO2 * X5 


ON 
Qu 
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